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ABSTRACT 

A disturbance observer achieves on-line estimation of unknown disturbances acting on a dynamic system based on 
the input and output information. This paper proposes a nonlinear disturbance observer for a class of nonlinear systems. 
Arbitrarily close estimate of the unknown disturbance can be obtained if the disturbance is bounded with bounded 
first-order time derivative, and the system nonlinearity satisfies a Lipschitz condition. A unique feature of the proposed 
nonlinear disturbance observer is that it allows the Lipschitz constant to be arbitrarily large. 
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1. INTRODUCTION 

In many engineering and science problems, one is concerned with the estimation of disturbances (or unknown 
inputs) acting on a dynamic system [1], [2], [3]. In control engineering, disturbance observer and unknown input observer 
are the most popular approaches to on-line disturbance estimation. The Disturbance Observer (DOB) [4] and [5] approach 
refers to using direct system inverse to estimate the input disturbance. In order to avoid a non-proper system inverse, 
a low-pass Q-filter is cascaded with the system inverse, where the Q-filter is chosen to achieve accurate disturbance 
estimation and to avoid amplification of measurement noise. Different design considerations and design skills have been 
proposed in the literature [6] - [11]. Disturbance rejection control based on DOB has been proved to be very effective 
especially in motion control applications. 

The second approach to disturbance estimation is the unknown input observer. In recent literatures, the most 
popular unknown input observer is the joint state and disturbance observer, where an observer is constructed to 
simultaneously estimate both the unknown state and the unknown disturbance. The designs in [12] - [14] are applicable 
under the assumption that the unknown disturbance is constant or slowly varying. The extended state observer design 
[15] - [17] basically follows the same assumption If knowledge of a disturbance model is known a priori, Luenberger type 
jointed state and disturbance observer can be constructed in [18] - [21] When no disturbance model is available, 
the designs in [22], [23] are applicable. 

Studies on the disturbance observer designs in linear systems are more common than nonlinear ones in the 
literature. The work in [22] and [24] simultaneous estimate the unknown disturbance and unknown state for a class of 
nonlinear systems. Their solution relies on solving an Linear Matrix Inequality (LMI), but the LMI admits a positive 
definite solution only if the Lipschitz constant of the nonlinearity is small [24]. Further, they constraint the relative degree 
from the nonlinearity to the system output to be one [22] There has been study of the disturbance observer for a particular 
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class of nonlinear system, that is, the robot systems [25] - [27]. In these literatures, the robot disturbance observer must 
assume that all the state variables (angles and velocities of all joints) are available for measurement. While the disturbance 
observer in this paper requires only information of system outputs (for example, angles of certain joints), and this 
requirement is much more relaxed than the robot disturbance observer designs. 

The disturbance observer design problem presented in the paper is more difficult than the state observer design 
problem for nonlinear systems. For the state estimation problem for nonlinear systems, one may see the interesting work in 
[28], [29], which shows that a nonlinear state observer exists as long as the distance to unobservability is larger than the 
Lipschitz constant. For the disturbance observer problem treated in this paper, one must first design a robust nonlinear 
observer that ensure small estimation error in front of unknown disturbances (Theorem 3 below). Then, one needs to prove 
that in the state estimation error dynamics, small state estimation error implies small input ex citation. After establishing 
this property can one prove that the disturbance estimation error is small. The disturbance estimation problem is therefore 
non trivial, and is dealt with by Lemma 4 and Theorem 5 in this paper. 

The remainder of this paper is arranged as follows. Section 2 presents the new nonlinear disturbance observer 
design for systems whose nonlinearity enters the system equation in the same direction as the disturbance. Section 3 further 
extends the nonlinear observer so that the nonlinearity and disturbance enter the system equation in different directions. 
Finally, Section 4 gives the conclusions. 

2. NONLINEAR DISTURBANCE OBSERVER 

Consider the disturbance estimation problem for a nonlinear system 

x - Ax + Bu + G(f(x) + d), ^ 

y = Cx. 

Where re R" is the unknown system state, u e R m the control input, ye R p the measured system output, 
d e R q the unknown disturbance satisfying 

\\d\\<D 0 , ( 2 ) 

For two positive constants n and l) l , and f(x) e R q is a nonlinearity satisfying the Lipschitz condition: 

| fix) ~ /(x)|| < y\\x - x\\, (3) 

For some Lipschitz constant y > 0 . Note thaFt in this paper the Lipschitz constant y can be arbitrarily large 
while previous papers all constrain y to be small. C and G are full rank. It is assumed that 

dim( y ) > dim( d ) , (4) 

Where dim stands for dimension. Without loss of generality it is assumed in the sequel that dim(y) =dim (d); the 
case dim(y)>dim (d) can be treated by artificially expanding the dimension of d. Further, it is assumed that (A, C) is 
observable, (A, G) controllable, and the square system (A, G, C) has only stable zeros (minimum-phase). 
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To solve the disturbance estimation problem for the nonlinear system (1), one constructs a nonlinear joint state 


and disturbance observer as follows. 

x — Ax + Bu + L(y - Cx) + Gf(x), (5) 

d = 4ftW{y-Cx), ( 6 ) 

Where d is a disturbance estimate for the unknown disturbance d, W e R pxp is an unitary matrix to be 
specified, the observer feedback gain L e R" xp is designed as in the LTR observer [30] and [31], 

L = QC T , (7) 

Q{A + aI)‘ +(A + aI)Q-QC T CQ + KGG T = 0 , 7r>0, a>0. ( 8 ) 


In the above equation, it is assumed that (A, G, C ) is minimum-phase; that is, all system zeros of (A, G, C) are 
in the open-loop left half plane. The design parameter a is chosen to be smaller than the minimum of | Re( z, ) | so that the 

assumption (A, G, C ) is minimum-phase ensures that (A + al, G, C ) is also minimum-phase. Previous studies [24] on the 
nonlinear disturbance observer require the solution of an LMI problem, which is equivalent to an type Riccati 

equation. However, their Riccati equation admits a positive definite solution matrix only if the Lipschitz constant y in (3) is 
small. In contrast, this paper employs an LQ type Riccati equation (8), which is y independent. The existence and positive 
definiteness of the solution matrix Q of the LQ Riccati equation (8) is always guaranteed [32] and [33] for all design 
parameters n > 0 and a > 0 if (A, C) is observable and (A, G) controllable. 

Remark 1 

Note that the above proposed robust observer is different from the linear high gain observer [34] because the 
proposed observer is nonlinear due to a nonlinear term f (x) in (5). Furthermore, even if a linear high gain observer can 
work in state estimation, there is no disturbance estimation mechanism in a linear high gain observer. 

The solution matrix Q £ R" x " in (8) depends on its design parameter n, and Q(n) and ;t satisfies the following 
relationship. 

Lemma 1 [30], [35] 

If the system (A, G, C) is square and minimum-phase, the solution Q(n) of the observer Riccati equation (8), 

satisfies Jim l - o. 

*-»“ n 

With the above lemma, one can further derive the following result. 

Lemma 2 

If the system (A, G, C) is square and minimum-phase, the observer feedback gain L in (7) satisfies 

-> GW as n -» °o, ( 9 ) 

dn 
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where W e R pxp is an unitary matrix satisfying WW 7 = I . 

Proof 

Divide the Riccati equation (8) by n to obtain 

^-(A + aI) T + (A + aI)^--—^-— + GG T =0. 
x x x 1 jn x 

It follows from Lemma 1 that the first and second terms in the above equation can be neglected as n — » oo, and this 

results in -> GG T . 

six s]X 

Solving this equation yields QC T -A sfxGW for some unitary matrix W. Finally, quoting (7) gives 
L — > sfxGW ■ End of proof. 

Remark 2 

The unitary matrix W in Lemma 2 can be obtained from (9) as follows. \y _> -J (G' G) ' G' L as X — > °o. 

six 

The state estimation error X = X — X resulting from the above proposed nonlinear observer satisfies 

x = {A-LC)x + G(f(x)-f(x) + d). (10) 

The goal of the following theorem is to prove that even for large Lipschitz constant in (3), the state estimation 
error X decays to almost zero in the error dynamics (10) if the design parameter n in (8) is chosen sufficiently large. 

Theorem 3 

Consider the nonlinear system (1) with (A, G, C) square and minimum-phase. For whatever large Lipschitz 
constant y in (3), the state estimation error resulting from the nonlinear observer (5) asymptotically satisfies 

lim llxll < £ , (11) 

oo 11 11 

Where £ can be an arbitrarily small positive constant if the design parameter n > 0 in the Riccati equation (8) is 
chosen sufficiently large. 

Proof 

Define a Lyapunov function V = x' Q X , where Q > 0 is from the observer Riccati equation (8), and x as in 
(10). The time change rate of V along the trajectory (10) satisfies 

V<-2cW-\^-x\G T Q- l x^+^f(x)-f(^+\d\)\G T Q-A\<-2cW-\yf-x\G T Q^+?lYl^+D 0 )\G r Q- l x\, 

Where y = Cx and the Lipschitz condition (3) was used to obtain the second inequality. Note that the maximum 
of the last two terms in the last equation occurs when ||G r <2“ 1 xj| = (y||jc|| + D 0 )/x , with the maximum value being 
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(y||x|| + Z) 0 ) 2 /;r • Hence, 


, ,,2 (r\\x\\+D„f (r|x| + Z> 0 ) 2 

V<-2cW-\y\ +HJLJ < — ctrK — («V" — °L_) 

71 71 


( 12 ) 


From the last inequality, one observes that V decays exponentially as long as V > (?'||jf|| + D 0 ) 2 j{0CK) ■ Therefore, 
eventually we have V < (y|ic|| + D 0 ) 2 /(a7r) as t -> <*>. Using V > a(|2 _1 )||ic|| 2 = |jf|| 2 /a(0, we obtain 


ic < 


cr(Q/ff) 


(y||jc|| + D 0 ) as f-> oo. Re-arranging the equation gives lim(l - y, 


g(Q/7T) 


)| x ( f )|<£ > 0 


cr(Q/ff) 


a V a V a 

Lemma 1 ensures that given any larege Lipschitz constant y, there always exists a sufficiently large observer design 
parameter n such that the number in the parenthesis in the above equation is positive. Hence, we can write 


D r 


lim ||x(f)|| < £ - 


a 


1 -Y 


a 


The above inequality together with Lemma 1 imply that given any bounded Lipschitz constant y and disturbance 
upper bound D 0 , I lx I will eventually become arbitrarily small as long as n is sufficiently large. End of proof. 


Compared with the nonlinear state observer in [28], [29], the nonlinear observer in Theorem 3 is robust against 
unknown disturbance; it achieves arbitrarily small state estimation error in the fact of large unknown disturbance. 
Furthermore, even facing with arbitrarily large Lipschitz constant, the proposed observer is still stable. 

In the formulation of this paper, the disturbance d enters the system through a constant vector channel G in the 
system equation (1). For certain nonlinear systems, for instance, the robot systems or other nonlinear systems in [29], 
the system has the structure that the input channel vector G(x) is state dependent. In this case, we can use state 
transformation so that in the transformed state, the new system has a constant disturbance channel vector. For example, 
a nonlinear system 

x = Ax + G(x)d, G(x ) = G x n x (x) + G 2 n 2 (x), 

Where G and G 2 are constant vectors, and n \ (x) and n 2 (x) are nonlinear functions of x. If we defines extended 
state z = [x, x n+i ], where x n+i = d, then the above system can be written as 

x- Ax + G 1 n 1 (x)x„ +1 + G 2 n 2 (x)x n+1 , (13) 

x n+l -d l , where d { = d. (14) 

The above new system can be viewed as a nonlinear system with nonlinearity «i(x)x n+ i and n 2 (x)x n+ i and an 
disturbance d\ with constant disturbance channel vector. We can then follow the procedure in the next section to design 
state observer and disturbance observer. 

The state estimation error dynamics resulting from the proposed nonlinear observer can be written as, via the 
result of Lemma 2, 
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x = Ax + G(f(x) — f(x) + d) — LCx = Ax + Gv, as 7T — (15) 
v = f(x ) - /(.£) + d - \[?rWCx - f(x) - f{x) + d —d. (16) 

We can view the above error dynamics (15) as a linear system subject to an input excitation v. According to 
Theorem 3, the system state x becomes arbitrarily small if the design parameter n of the nonlinear observer is chosen 
sufficiently large. One question stemming from equation (15) is that since the state x is small, can one conclude that its 
input excitation v is also small ? This is the bounded state bounded input (BSBI) property. A counter example of BSBI 
property is the sliding mode control applied to a linear system, x = Ax + B{u + d), U : = —p ■ sign(s), where the 
system is subject to an unknown disturbance d, p is the switching control gain, and is the sliding variable. If ,v and p are 
properly chosen, the switching control u can drive the system state x to zero. In this case, the state x is small, but the input 
excitation u + d is not necessarily small. 

To answer the question of BSBI property, one quotes a well-known lemma in the study of Adaptive Control. 
Define the truncated norm of a time signal v as ||v|| = SUpj,< T < ( ||v( y r)|| . A time signal v is regular if 

||v|| < Mj ||v|| + M 0 for two positive constants Mj and M Q for t > T. Due to (11) in Theorem 3, there exists a finite 
time T such that 

||Jc(f)|| < 1.5f for all t>T. (17) 

In the sequel, this T will be used in the definition of ||v|| = SUp r < T < f ||v(T )|| . Hence, it follows from (17) 

that 

||jc(f)||<1.5e as r— > oo. (18) 

Note that the disturbance observer design problem for nonlinear systems is more difficult than the state observer 
design problem for nonlinear systems. In the disturbance observer problem, one must first design a robust nonlinear 
observer that ensure small estimation error in front of unknown disturbances. Then, we need to prove for the state 
estimation error dynamics that small state estimation error implies small input excitation. After establishing this property 
we can prove that the disturbance estimation error is small. The latter job is non-trivial, and is dealt with by Lemma 4 and 
Theorem 5 below. 

Lemma 4 [36] 

Given a linear system x = Ax + Gv , where (A, G ) is controllable, if v is regular, then 

||v||,< tf||jc||,, (19) 

For some positive constant K 

Lemma 4 states that as long as the input is smooth in the sense of being regular, we can use (19) to conclude BSBI 
property. With Lemma 4, one is now in a position to prove the main result of this section that d in (6) is a close estimate 
of the unknown disturbance d. 
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Theorem 5 

Consider the nonlinear system (1) subject to the unknown disturbance d , and the corresponding nonlinear observer 

(5) and the disturbance estimate d in (6). If functions r H— ^JL /\ x q f ( x ) , are all Lipschitz, and the control input u 

dx dx dx 

is bounded, the disturbance estimate d in (6) is arbitrarily close to the true disturbance d if the observer design parameter n 
in (8) is sufficiently large. 

Proof 

Define a function 

M(x,t) = —Ax + — Bu + — G(f(x) + d). 
ox ox dx 

From hypothesis of this theorem, M is Lipschitz in x. According to the definition of 

v = /(x) - f{x) + d - yfftWCx in (16), the time derivative of v can be expressed as 

v = M(x,t)-M(x,t) + d- 4 ttWC{Ax + Gv ). 


Finally, quoting the Lipschitzness of M(x,t), boundedness of d , and Theorem 3, one concludes from the above 
equation that v is regular. Since v in (15) is regular, we deduces from Lemma 4 that 
||v(f)||<||v|| ( < X'||jc|| f <1.5A"f as t — > oo_ where the last inequality follows from (18). Define the disturbance estimation 

error d = d — d , and re-arrange (16) to yield 


\d 


d-d 


|v — /(x) + /(x)|| < ||v|| + )'||x|| < 1 .5 Ks + yE, where we have used the Lipschitz condition (3) to 


obtain the first inequality. Since e approaches zero as the design parameter n approaches infinity, the last inequality shows 
that the disturbance estimation error d can be made arbitrarily small. End of proof. 


The above theorem shows that the disturbance estimate d is arbitrarily close to the true disturbance d if the 
observer design parameter n is sufficiently large. In real applications, we can use computer simulations to test just how 
large the design parameter n should be so that a desired disturbance estimation error is obtained. 


3. EXTENSION TO MORE GENERAL SYSTEMS 


In previous sections, the unknown disturbance and the nonlinear function enter the state equation through the 
same channel. The purpose of this section is to construct a disturbance observer for systems whose unknown disturbance 
and nonlinear function enter the state equation through different channels. Consider a nonlinear system subject to an 
unknown disturbance, 

x = Ax + Bu + GJ(x) + G 2 d, 

( 20 ) 

y = Cx. 

Where all symbols are as defined in (1), and G, G, ■ Further, it is assumed that (A, C ) is observable, and 
(A, [Gi, G 2 ]) controllable, f (x) G R ql , d e R ql , and y G R p . 
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Constructing a disturbance observer for the above system is more difficult than for the system (1) in the previous 
section. However, it will be shown that if the number of output sensors is sufficient in the sense that 

dim(y ) > dim(f(x )) + dim(d) O p>ql + q2 , (21) 

Where dim stands for dimension, then design of a nonlinear disturbance observer is achievable. For simplicity, it 
will be assumed in the sequel that dimiy) = dim{f{x )) + dim(d) , rank(C) = rank([G 1 ,G 2 ]) . 

The case dimiy) > dim if (x)) + dim(d) can be treated similarly by artificially expanding the dimension of d 
or fix). Note that this condition (21) is similar to the assumption adopted in other reference [24]. 

A nonlinear disturbance observer for the system (20) is proposed as follows, 

x - Ax + Bu + L(y - Cic) + G { f (x), (22) 

d = [0 ?2X?1 , I q2xq2 ]yfzw ( y-Cx ), (23) 


Where W is an unitary matrix to be specified, and the observer feedback gain L e R mp is given by 

l = qc t , 

Q(A + aI) T +(A + aI)Q-QC T CQ + x[G l ,G 2 ][G l ,G 2 f = 0 , x>0, a> 0 . 


(24) 


Note that the above Riccati equation is different from that in the previous section. The existence of a positive 
definite solution matrix Q G R" x " is guaranteed [32] and [33] when (A, C) is observable, and (A, [G \ , G 2 ]) controllable. 

The state estimation error x = X — X resulting from the proposed robust nonlinear observer (22) obeys the 

/(*) - /(*) 


dynamics x = (A - LC)x + [G l , G 2 ] 


where the last term satisfies an upper bound 


f{x)~ f{x) 
d 


< + q < y % + D (t in which we have used (3) and (2) to obtain the second inequality. 


The following theorem shows that the above robust nonlinear observer (22) guarantees that the estimation error 
X — X decays to almost zero even in the face of large disturbance d. 

Theorem 6 

Consider the nonlinear system (20) subject to a bounded disturbance d , and the assumption that (A, [Gi , G 2 ], C) 
is square and minimum-phase. For whatever large Lipschitz constant y, and whatever large disturbance d, the robust 
nonlinear observer (22) guarantees that the state estimation error x — x—x will eventually approach a small residual set 
containing the origin, with the size of the residual set approaching zero as the design parameter n in (24) approaches 
infinity. 

Proof 

The proof is similar to that of Theorem 3, and is omitted. 
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The next theorem proves that the disturbance estimate d in (23) is arbitrarily close to the true disturbance d if the 
design parameter n in the Riccati equation (24) is chosen sufficiently large. 

Theorem 7 

Consider the nonlinear system (20) subject to the unknown disturbance d, and the corresponding nonlinear 

observer (22). If functions Ax, Gf{x )> are all Lipschitz, and the control input u is bounded, the 

dx dx dx 

disturbance estimate d in (23) can be arbitrarily close to the true disturbance d if the observer design parameter n in (24) is 
sufficiently large. 

Proof 


The proof is similar to that of Theorem 5, and is omitted. 

To demonstrate the effectiveness of the proposed nonlinear disturbance observer, we present a numerical 
simulation example. 

Example 



~ 0 

1 

0 

o' 


'O' 


'O' 


'0 

f 


-1 

0 

1 

0 


0 


1 


0 

0 

A = 





, B — G 2 — 


, G,= 


, c r = 




0 

0 

0 

1 


0 


0 


1 

0 


1 

0 

1 

-1 


1 


0 


0 

0 


Consider a single-link with flexible joint robot [37], whose governing equation is as in (20) with the control u(t) = 
0, the disturbance d{t) = 2, and nonlinearity fix) = — I Osin A' . Note that the nonlinearity fix) is globally Lipschitz, and the 

large Lipschitz constant y = 10 implies a long or heavy robot arm. The initial condition are x(0) = [1, 1, 1, 1] for the 

robot, and X (0) = [-1, -1, -1, -1] for the observer (22). The modified Riccati equation (24) has design parameters 

CC — 1 , and ir(t) scheduled continuously from 1 to 10 12 within 10 seconds according to the formula n (t) = io 12t/10 , t £ [0, 
10] and n{t) = 10 12 , t > 10. Note that the design parameter k is time-varying only for the initial transient 10 seconds. After t 
>10 seconds, the design parameter n becomes constant. Hence, the analysis in the previous sections applies Theorems 3 
and 5 holds for t > 10. At the steady state, n = 10 12 , the corresponding observer feedback gain is found to be 


L = 


7.071 le-004 1.4162e + 003 
1.0004 1.0014e + 006 

1.4152e + 003 7.071 le- 004 
1.0000e + 006 1.0004 


Figure 1 shows the disturbance d and the proposed disturbance estimate d . It is seen that the d converges to d, 
as predicted by Theorem 5, in about 4 seconds. Figure 2 shows the state estimation error ||x(/)|| versus time. 

The estimation error decays to almost zero (||x(f)|| = 10 3 ), confirming that the proposed robust nonlinear observer 
design is successful in the face of disturbances and large-Lipschitz-constant nonlinearity. Note that if we use n = 10 12 right 
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from the beginning of observation, there will be a substantial peaking phenomenon. The proposed stepwise scheduling of n 
(t) has successfully avoided the peaking phenomenon. 



Figure 1: Disturbance and Its Estimate 



Time (sec) 

Figure 2: Norm of State Estimation Norm 

4. CONCLUSIONS 

Previous studies on the disturbance observer are mostly linear systems. This paper proposes a nonlinear 
disturbance observer for a class of nonlinear systems. Arbitrarily close estimate of the unknown disturbance can be 
obtained if the disturbance is bounded with bounded first-order time derivative, and the system nonlinearity satisfies a 
Lipschitz condition. The contribution of this paper is that it allows the Lipschitz constant to be arbitrarily large, and it 
imposes no constraint on the system relative degree. 
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